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Spacecraft historically have had sub-lkW e electrical requirements for GN&C, science, 
and communications: Galileo at 600W e , and Cassini at 900W e , for example. Because most 
missions have had the same order of magnitude power requirements, the Power Distribution 
Systems (PDS) use existing, space-qualified technology and are DC. As science payload and 
mission duration requirements increase, however, the required electrical power increases. 
Subsequently, this requires a change from a passive energy conversion (solar arrays and 
batteries) to dynamic (alternator, solar dynamic, etc.), because dynamic conversion has 
higher thermal and conversion efficiencies, has higher power densities, and scales more 
readily to higher power levels. Furthermore, increased power requirements and physical 
distribution lengths are best served with high-voltage, multi-phase AC to maintain 
distribution efficiency and minimize voltage drops. The generated AC-voltage must be 
stepped-up (or down) to interface with various subsystems or electrical hardware. Part of 
the trade-space design for AC distribution systems is volume and mass estimation of high- 
power transformers. The volume and mass are functions of the power rating, operating 
frequency, the ambient and allowable temperature rise, the types and amount of heat 
transfer available, the core material and shape, the required flux density in a core, the 
maximum current density, etc. McLyman [1] has tabulated the performance of a number of 
transformers cores and derived a “cookbook” methodology to determine the volume of 
transformers, whereas Schawrze [2] had derived an empirical method to estimate the mass 
of single-phase transformers. Based on the work of McLyman and Schwarze, it is the intent 
herein to derive an empirical solution to the volume and mass estimation of three-phase, 
laminated El-core power transformers, having radiated and conducted heat transfer 
mechanisms available. Estimation of the mounting hardware, connectors, etc. is not 
included. 
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Nomenclature 

core width [m] 

ratio of window width, window height, and core depth, respectively, to core leg width 
core cross-sectional and window areas [m 2 ] 

Area product [m 4 ] 

flux density [Tesla] 

flux [Weber] = <f> sin(cot) 

permeability of free space, 4jt x 10" 7 [T-m / A] 

radiating surface area [m 2 ] 

shielding (view) factor 

Stefan-Boltzman constant, 5.67 x 10" 8 [W / (m 2 -K 4 )] 
emissivity of surface 

transformer maximum, ambient, and cold-plate temperatures, respectively [K] 


I. Introduction 

P OWER transformers are used to change distribution voltage levels or to electrically isolate one system from 
another. The size and mass of the transformer is a function of many design trade-offs (e.g., core 
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permeability/size versus efficiency) and design points (e.g., maximum operating temperature, power rating, and 
frequency). This paper is not intended to be a thorough transformer design treatise, but rather it presents an 
empirical solution to estimate the volume and mass of a three-phase, El power transformer, which has only radiation 
and conduction heat transfer available and will operate at no more than 10kHz. The method can be extended for 
other core geometries, material, as frequencies. 


T 


II. Transformer Equations 

o estimate the volume of a transformer, the electrical, 
mechanical, and thermal performance for a given set of 
criteria must be derived. 


A. Electrical Equations 

The general voltage equations for a gapped, laminated EI- 
core, for which the dimensions are given in Figure 1, are derived 
for each phase. A voltage on the primary coils of each phase 
induces a time-varying flux, ^(t), into the core. This flux links 
with the other windings and induces voltages in the secondary 
windings. Each phase winding is on an individual leg (i.e., one 
primary and one secondary per leg) and equally shares the 
window area with another phase. Phases 2 and 3 lag phase 1 by 
120° and 240°, respectively. With the gaps being equal, the 

reluctance of each leg is equal and, therefore, the one-half of the flux induced by a primary winding 
other two legs. The resultant voltage equations for any winding of n-turns are given in Eq. (1). 



Figure 1. Mechanical dimensions and flux flow 

flows in the 
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Simplification of Eqs. (1) by use of the following identities 

</>(t) = <&sm(cot) 
d 
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cos(a-b) = cos(a)cos(b) + sin(a)sin(b) 
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results in expressions for the voltages: 


V' 1 (f) = 1.5« 1 <t>(D cos (cot) 
v 2 (t)=l.5n 2 <t>co cos (or-—) 


Art 

v 3 (t) = 1.5 n 3 T> co cos( co t ) 


(2a) 

(2b) 

(2c) 


From Fig. 1 it is seen that the cross-sectional area, A c , of the core equals da 2 . For a given core material the 
maximum, non-saturating flux density, B , should not be exceeded. By definition <j> = b * ,4 . and the root-mean 
square (rms) value of the primary and secondary voltages, using Eq. (2), can thus be expressed as 


1.5 N p O)BA c 
1 .5N s O)BA e 


(3a) 

(3b) 
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B. Fundamental parameter derivation 

The fundamental dimension, a [m], of the transformer core must be found. Equation 11.10 in McLyman’s handbook 
relates the current density of the transformer to the area product. No distinction is made between primary or 
secondary current. Thus the rms-current densities are implied to be equal, J = =J [A/m 2 ]. This implies that 

each window area, W a , must be shared equally by the windings within it. Because the primary and secondary of two 

W 2 

phases share the window, we have: w = w = — [nr]. 

A r a,p f a, s ^ 

The ideal cross-sectional area of a winding is a function of the number of turns, N, and the current density 
allowable within the window area (all four equal). Only a portion of the window can be filled with copper because 
of insulation and inter-winding gaps. This reduction is denoted as the window utilization factor, k u . Thus, the ideal 
cross-sectional area of a single turn of a primary and a secondary winding is: 


A l = k ' ,Wa 
4 N p 

[m 2 ] 

(4a) 

A - KW ° 

- [m 2 ] 

(4b) 

2 4 N, 




Because Eq. (4) will probably not result in an actual AWG wire areas, the transformer design must be adjusted to 
accommodate the next-larger, actual AWG diameter, based on the wire diameter relationship for gage sizes stated by 
Fink & Carroll [3], 

The primary and secondary currents are fixed by the phase voltage and by the rated power of the transformer. A 
further constraint, which must be met, is based upon the relationship between J and the cross-sectional area of the 
conductor. Namely, I = J A conductor , which leads to: 

, Jk„W a rA n 

/ = [Arms] (5a) 

" 4 N P 

I [A rms ] 

1 4 N s (5b) 

The power of a three-phase transformer, operating in a balanced mode, is simply three times the power in any 
phase (i.e., S rms = 3 I rms V m s). By substituting the voltage equations, Eq. (3), the current equations, Eq. (5), and 
McLyman’s definition for the area product (i.e., A p = 1.5 W a A c [m 4 ]) into the power equation results in: 

„ „ Jk„W„ 1.5 NaBA c 0.75 , , , - fVAl 


S = 31 V =3 : 

P P P 


-j^-A (JaBk u ) [VA] 
V2 


S =31 S V S =3 J k “ Wa L5N ’“ BA ‘ =™A,iJ0*k.) [VA] 

5 s 4N S V2 ( bb ) 

For an ideal, loss-less transformer, the output power equals the input power and, thus, either expression in Eq. 
(6) can be used to solve for A p . The resulting expression for A p as a function of the rated power, electrical 
frequency, actual current density within the conductors, the core utilization factor, and the peak flux density is: 


A 

3 JcoBk 


m,ed [m 4 ] 


Having found an expression for A p , the fundamental core dimension “a” can be found as follows: 

A p =L5W a A c [m 4 ] 

= 1.5 (hk a 2 ) (d a 2 ) 


A p = l.5d hka 4 [m 4 ] 
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C. Copper and core losses 

The transformer losses are the sum of the copper and core losses. First, the copper losses will be determined. 
By definition R = pL/A, where A is the cross-sectional area [m 2 ] of a conductor, p is the conductor’s electrical 
resistivity [fi-m], and L is the total length [m]. Having found the cross-sectional area of the winding per Eq. (4), the 
total length of a winding must be found. To accomplish this, the Mean Length of a Turn, MLT, of an n-turn winding 
(whether primary or secondary) is found according to: 

MLT = Total length = l max+ l mm [m] 

turns 2 ' ’ 

From Fig. 2 the MLT can be approximated as follows. Assume that a O-diameter wire is bent at 90° angles, then the 
innermost [] min ] and outermost (l max ) winding lengths are: 


l mm =2c,a + d) [m] 

l mm =2a{(k+d) + (k + \)} = 2a{2k+d + \) [m] 

Substituting these lengths and Eq. (8) into Eq. (9) and simplifying results in an expression for the MLT: 

MLT = 2a(k + d + 1) 

= 2(k + d + l)4p 51 [m] (10) 
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Figure 2. Transformer dimensions with windings. 


Now the total copper loss, ignoring skin effect in this derivation, for the three phases is: 


P Cl =3(I;R p +I;R s )[W] 


= 3 
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(11) 


The core loss is the sum of the hysteresis and eddy current losses, and can be approximated over the range of 
frequencies of interest as [4]: 
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Paresis =Volume core k h B' 6 f [W] 

P edd y=V°lume core k e t 2 f 2 B 2 [W] 

The lamination thickness, t [m], is set at 90% of the skin depth, which is a function of frequency and of the 
relative permeability and resistivity of the material. The minimum value allowed is 0.007” (0.18 [mm]), and is found 
according to Eq. (14). 


( 12 ) 

(13) 


(14) 


*(/) = >0.18 [mm] 

\XfMrM 0 

The core volume is found from Fig. (2) to be: 

V core = (ha + 2a)(3a + 2ka)d a—2(k a)(h a)(d a) = d ( 3h+4k + 6)a 3 [m 3 ] 

Substitution for the core volume into Eq. (12, 13), summing, substituting Eq. 8, and simplifying leads to the core 
power loss to be: 


Pcore =d(3h+4k + 6) 


r 2A „ ) 

|4 ( 

p 

f 

v 3i ihk j 

V 


k h B L6 +k e t 2 f B 2 


[W] 


(15) 


The electrical equations are complete. However, if these were used to solve the transformer volume, without 
considering the heat dissipation and maximum temperature limitations of a physical device, the solution would be an 
unrealistically small transformer. Therefore, the heat rejection must be estimated. 

D. Heat Rejection Equations 

A transformer's size is a function of its efficiency (power loss) and the ability to reject heat to the environment, 
among other variables. Heat transfer methods are radiation, convection (moving fluids), and conduction (solids). 
For space applications convective heat transfer to the surroundings is non-existent. For the volume estimation 
problem, heat rejection is thus limited to radiation and conduction through the core to the cold plate. Furthermore, 
the surface temperature of the cold plate is assumed to be constant. 

The heat transfer for non-black body surfaces is given by: 

P md =F [W] (16) 

It is the assumed that the transformer is enclosed with a shielding factor of 0.9. Transformer surfaces are 
assumed to not radiate heat to other surfaces of the transformer. The exposed area of the transformer is defined as 
the total surface area less the surface area in contact with the cold plate. From Fig. (2) the total surface area of the 
transformer can be found: 

A S a = 2(ka + da)(3ka + 3a) + 4(2ka + 3a)a + 4da" + 2ha(3ka + 3a) + 2ha(ka + da) 

A<. a = 2 ct (3k 4- Ik + 3dk + 3d -t- 6 -t- 4-hk -t- 3h -t- hd ) [m~] 

And the area of the core in contact with the cold plate is 

Ao P _core =d a(2k a + 3a) = a 2 (2d k +3d) [m 2 ] 

By subtracting Eq. (18) from Eq. (17), the surface area available for radiation heat transfer is 

A r = a 2 (6k 2 + 14k + 4dk +ld + 1 2 + 8 hk + 6h + 2 hd) 

A r = p 4 ^- (6k 2 + 1 4k + 4dk +ld + 1 2 + 8 hk + 6 h + 2 hd ) [m 2 ] 

V 3dhk 

As a result, the expression for radiated heat transfer becomes: 


(17) 

(18) 

(19) 
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( 20 ) 


p md = F (6k 2 +14k+4dk+ld + l2+8hk + 6h + 2hd)£(r{T^ - T a 4 ) [W] 

V 3<7/tk 

Because the transformer surfaces will not have convective heat 
transfer available (i.e., no moving air), the cold plate must transfer 
the heat flux, which is not radiated. Cold plates have two heat 
transfer modes: conduction and convection. Heat is conducted 
through the core, across interface, and through the cold plate to the 
coolant interface. Heat transfer from the cold plate material to the 
moving coolant is via convection. This is not a cold plate design. 

Therefore, it will be assumed that the cold plate is designed to 
maintain the cold plate surface at a constant and uniform 
temperature, T cp . To further simplify the problem, it will be 
assumed that all of the heat flux conducted within the transformer 
begins at a surface area of uniform temperature (the maximum 
allowed). This heat flux is uniformly distributed at one end of an equivalent solid core, having no windows areas 
and has a height equal to the mean height of the core as shown in Fig. (3). All heat flux must pass through this 
volume to the cold plate. Lastly, the contact surface area is that given in Eq. (18). 

The conducted heat transfer capacity, P con d, of this equivalent core mass with a given thermal conductivity, x c 
[W/m-C], is given by: 

p _ ^max ~ T cp ) A cp (21) 

mean height 



Figure 3. Equivalent core mass for heat 
transfer 


Substituting for the bottom surface area of the cold plate and for the mean height, the expression for the thermal 
energy transferred by the core to the cold plate is: 


5 (^max T C p ) ^4 


P =- 

cond 


2 Ap 

3 dhk 


*(2dk+3d) 


3/t+10 


[W] 


( 22 ) 


The thermal losses must equal the heat transfer system design capacity to limit the fully loaded transformer to its 
maximum operating temperature. This can be written as 

Pcore Pcu — Pcond Prad [W] (23) 


E. The system of Equations 

The transformer problem consists of a system of equations (i.e., Eqs. 7, 11, 15, 20, and 22) subject to further 
constraints (i.e., Eq. 23, h > 0, d > 0, k > 0, ./ < J max ). From the solution of this system of equations all other 

parameters can be found. For example, “a” is found according to Eq. (8), MLT by Eq. (10), Np from Eq. (3a), etc. 
Also, having found N p and N s , the wire gage can be determined according to Eq. (3a) and the transformer ratio, since 
2 primary and 2 secondary windings share the window and have a given utilization factor. Now all of the 
parameters are available to estimate the mass, dimensions, and volume of the transformer copper and core as a 
function of it power, voltage, and temperature constraints. 


III. Sample Estimation Results 


A 50kVA, three-phase, 400V:400V 
isolating transformer is designed to operated with 
1.5[T], J max = 3.5 x 10 6 [A/m 2 ], and the 
laminations are a function of frequency, but 
limited to a minimum of 0.18 [mm]. The core has 
parameters of p r = 1000, p = 9.579 x 10' 8 [fl-m], 
and a density of 7300 [kg/m 3 ]. The windings are 
copper for which p = 1.724x 10" 8 [fl-m] and 
density = 8890 [kg/m 3 ]. The window utilization 


Table 1. 50kVA transformer volume and mass Estimates 


f [Hz] 

50 

1000 

1500 

volume [cm A 3] 

34525 

9947 

8887 

Length [mm] 

416 

468 

458 

Height [mm] 

316 

192 

103 

Depth [mm] 

262 

110 

187 

mass [kg] 

169 

28 

25 
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is assumed to be 40%. The ambient, cold plate, and maximum- 
operating temperatures are 303K, 293K, and 333K, respectively. 

The transformer mass and dimensions — excluding mounting 
hardware, etc. — is shown as a function of frequency in Fig. (4). 
Furthermore, a comparison of the estimates for the 50kVA 
transformer operating at 50FIz, 1000Hz, and 1500Hz are given in 
Table 1. 


IV. Conclusion 

An empirical solution has been derived to estimate the volume 
and mass of cut-core, E-I, high-power transformers to be used in 
space applications, and have only radiation and cold-plate cooling 
available and must be sized to operate within the prescribed 
operating limits. Mounting hardware and connectors are not 
estimated, and the estimation is not an optimization of the 
transformer design. Furthermore, using the derivation steps in this 
paper, it should be possible to adapt the procedure to other core 
shapes and materials. 
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Figure 4. 50kVA transformer dimensions and 
volume 
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